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Abstract 

We show that in the late-infall model of our galactic halo by P. Sikivie 
the expected phase of the annual modulation of a WIMP halo signal in 
direct detection experiments is opposite to the one usually expected. If a 
non-virialized halo component due to the infall of (collisionless) dark mat- 
ter particles cannot be rejected, an annual modulation in a dark matter 
signal should be looked for by experimenters without fixing the phase a- 
priori. Moreover, WIMP streams coming to Earth from directions above 
and below the galactic plane should be expected, with a characteristic 
pattern of arrival directions. 



1 Introduction 

The dark halo of our galaxy may consist of WIMPs (weakly interacting massive 
particles). Direct detection experiments attempt to measure the nuclear recoil 
caused by these dark matter particles interacting with the material in a detector. 
The most important signature of a halo dark matter signal in these experiments 
is its annual modulation [pj. The WIMP interaction rate, the energy deposited 
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per collision, as well as the annual modulation of the signal depend strongly 
on the velocity distribution of the halo dark matter particles with respect to 
the detector. It is therefore very important to explore the possibility of halo 
particles velocity distributions which differ from the one usually assumed. The 
standard assumption is of a virialized dark halo, on average at rest in the rest 
frame of the galaxy, with a gaussian velocity distribution truncated due to the 
escape velocity from the galaxy, of about 600 km/s. The average WIMP velocity 
on Earth in then due to the motion of the Earth with respect to the galaxy (at 
about 200 km/s). This velocity is maximal around June 2 each year, when the 
velocity of the Earth around the Sun adds up maximally to the velocity of the 
Sun with respect to the galaxy, and minimal six months later, around December 
2 each year. 

A non-standard halo model is the late-infall model of Refs. 0, modified 
in recent years by P. Sikivie and collaborators 0, ||, 0|. It assumes a non- 
virialized dark halo in which collisionlcss dark matter particles falling into the 
galaxy oscillate in and out many times. At a given location in the galaxy, 
multiple flows of particles are possible, each having a specific velocity in a specific 
direction. The non-virialized flows of dark matter particles produce a velocity 
distribution completely different from the standard truncated Gaussian. The 
particular late-infall model of Sikivie in Ref. || is a self-similar axially symmetric 
infall model with net angular momentum and parameters adjusted to describe 
well our galaxy. For this model, which we call Sikivie's LI model or SLI model 
from now on, the local velocities and densities of the first twenty pairs of flows 
are given in Table 1 of Ref. We reproduce this table in our Table 1 for 
convenience. The first pair of flows corresponds to particles coming into the 
galaxy for the first time from opposite sides of it, the second to those passing 
for the second time, etc. 

Here we would like to clearly expose the differences that this particular late 
infall halo model by Sikivie implies for direct dark matter detection experiments. 

2 Characteristics of the WIMP Wind 

In the first comparison between the standard halo model and Sikivie's late infall 
model, we consider the flux of WIMPs as a function of solid angle for an observer 
moving with the Sun. If v is the velocity of a WIMP with respect to the galaxy, 
its velocity u with respect to the Sun is simply given by a galilean transformation 

u = v-v Q , (1) 

where 

v =v LS r + Vq. (2) 

v Q is the velocity of the Sun with respect to the galactic rest frame, vlsr is 
the velocity of the Local Standard of Rest (LSR), which is in the direction of 
galactic rotation, and v' is the peculiar velocity of the Sun, i.e. its velocity 
with respect to the LSR. In a coordinate system in which X points toward the 
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Table 1: Local densities and velocities of the first 20 pairs of flows in the self- 
similar infall model of ref . [|| . The first five columns are from Table 1 in Ref. @. 
The last column gives the flow speeds with respect to the Sun (as obtained by 
us). 
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galactic center, Y toward the direction of galactic rotation, and Z toward the 
North Galactic Pole, we adopt 

v LSR = (0, 220, 0) km/s, (3) 

for the velocity of the Local Standard of Rest, and 

v Q = (10, 13, 7) km/s (4) 

for the Sun peculiar velocity. (The uncertainty in the Sun peculiar velocity is 
of the order of 0.2 km/s in the Z direction and of as much as 3 km/s in the X 
and Y directions || ; the corresponding uncertainty in the phase constant of the 
modulation we discuss below is of several days.) 

With these preliminaries, the flux of WIMPs arriving at the Sun from within 
the solid angle dD, around the direction n is 

4c = Z- ( uf e (-un) ■ u 2 du. (5) 
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Here p is the local halo density, m is the WIMP mass, and fQ(u)d 3 u is the 
fraction of WIMPs with velocities with respect to the Sun within d 3 u around u. 
The WIMP velocity distribution in the rest frame of the Sun /©(u) is related 
to the WIMP velocity distribution in the galactic rest frame /(v) through 

/ (u) = /(u + v Q ). (6) 

Below we plot the WIMP flux d$ /dfl as a function of the arrival direction n 
in galactic coordinates (l,b), n = (cos b cos/, cos b sin/, sin 6). We do this for two 
models of the velocity distribution: the standard halo model and the late-infall 
model of Sikivie and collaborators. 

In the standard halo model, the WIMP velocity distribution is assumed to 
be a gaussian with velocity dispersion vq/V^ truncated at the escape velocity 

v2/w «, for|v|<U e 



/std(v)= 7V CSC ^3/ 2TJ 3 > -|v,^c S c (?) 

1 0, otherwise, 
where p is the local WIMP density, m is the WIMP mass, and iV csc = crf(zo) — 
2zq exp(— z^j/ir 1 / 2 with zq = v csc /vo is a normalization factor. For the sake of 
illustration, we take Uo = 220 km/s and v csc — 650 km/s. Other values do not 
change our conclusions. 

A simple integration gives the flux per unit solid angle as 

d ®std pv Q 

-dnT = — ho{1 ^ (8) 

where vq — |v Q | and 



y 2 o) 



-V^e y oy (3 + 2yl)[eri{^zl + yl-xl) - erf(y )] 
-2e^(l + 4y 2 + z 2 -x%~ 3y yfz$ + y 2 - x 2 )}. (9) 
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with x = vq/vo, 2/o = -v© • h/v , and z = v csc /v . 

In the late-infall halo model SLI, the WIMPs belonging to stream i have 
all the same velocity Vj and they contribute a density pi to the local WIMP 
density. The values of Vj and pi are given in Table 1. The model of Ref. Q 
does not include a velocity dispersion of the dark matter particles in the flows. 
In this case, the WIMP velocity distribution in the SLI model is 

/sli(v) = -£>£(v-v*). (10) 

P i 

The corresponding angular distribution of the WIMP flux simply consists of 
isolated points, 

d^shi sr^ Pi U i xi- l -M m\ 

i 

Here = Vj — v©, Ui = and = Ui/ui. If we introduce by hand a velocity 
dispersion Hi / for each flow, and assume a gaussian distribution of velocities 
for each flow, the velocity distribution becomes 

/sLTW=E^=3 e ^' (I 2 ) 

and the WIMP flux per unit solid angle 



— W,b), (13) 



where 

-x 2 

b) = 2(1 + yf) - v^Fe»? yi (3 + 2yf)[l - erf( W )]} (14) 

with Xi — Ui/vi and yi — Ui ■ n/vt. In the examples below, we set vi = 30 km/s. 

Figs. 1 and 2 show sky maps of the flux of dark matter particles in the 
standard halo model (Fig. 1) and in the SLI halo model (Fig. 2). These are the 
fluxes according to an observer at rest with the solar system. The sky maps 
in Figs. 1 and 2 are equal-area projections of the celestial sphere in galactic 
coordinates^] The galactic center is at the center, the galactic north (south) 
pole is at the top (bottom) , the galactic plane is horizontal passing through the 
center of the map. The Sun is moving towards the direction indicated by the 
white cross (Sun apex) . The white circle indicates the direction opposite to the 
Sun's motion (Sun anti-apex). The color in the maps shows the intensity of 
the WIMP flux coming from each direction (a lighter color indicates a larger 
intensity) . 

lr They are Mollweide projections. The relationship between the coordinates (x,y) on the 
map and the galactic coordinates (b, I) is x = -(2 v / 2«cos6»)/7T, y = \/2 sin0, where 9 is 
given by 29 + sin(2S) = 7rsinb. The inverse formulas are b = arcsin {[(29 + sin(20)]/7r}, 
/ = — (kx) / (2\/2 cos 9), where 9 = arcsin(y/v / 2). 
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Figure 1: Sky map in galactic coordinates of the WIMP flux as seen by an 
observer moving with the Sun for the standard halo model with v — 220 km/s 
and v esc = 650 km/s. Lighter colors indicate larger flux intensities. The WIMP 
wind comes mostly from the direction of the Sun motion (white cross) . 



Figure 2: Sky map in galactic coordinates of the WIMP flux as seen by an 
observer moving with the Sun for Sikivie's late-infall halo model, smoothing 
each stream with a gaussian with velocity dispersion Vi = 30 km/s. Lighter 
colors indicate larger flux intensities. WIMPs come from the direction of the 
streams (the bright spots). The most intense streams lie around the direction 
opposite to the Sun motion (white circle) . This reverses the phase of the annual 
modulation of a halo WIMP signal due to the motion of the Earth around the 
Sun. 
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In the standard halo model of Fig. 1, most WIMPs come from the direction 
towards which the Sun is moving. In the SLI model of Fig. 2, WIMPs come from 
the direction of the streams (the bright spots). In the figure, each stream has 
been convolved with a Maxwellian velocity distribution with velocity dispersion 
Ui = 30 km/s (otherwise the WIMPs of each stream would come only from 
one point in the sky). Fig. 2 clearly shows that in the SLI halo model most 
WIMPs come from directions in the hemisphere opposite to the Sun motion. 
As a consequence, the average 'WIMP wind" velocity on Earth, as pointed out 
by Sikivie himself in Ref. |6|, is reversed with respect to that in the standard 
halo model. This is because most halo particles in the SLI model move in the 
direction of the Sun's motion with a speed larger than that of the Sun. Therefore 
in the SLI model the annual modulation of a galactic WIMP signal has a phase 
opposite to the one usually assumed. Notice in Fig. 2 that in the SLI halo 
model many WIMPs come to Earth also from directions above and below the 
galactic plane. The directions of the streams form a characteristic "diamond" 
or "quad" pattern around the Sun anti-apex: the most intense WIMP fluxes 
are concentrated on the galactic plane and on the I — 270° meridian. 

3 Annual Modulation of WIMP Fluxes 

The time dependence of the galactic WIMP signal arises from the annual varia- 
tion of the Earth velocity with respect to the Sun. Direct detection experiments 
lacking directional capabilities are sensitive to the WIMP flux integrated over 
the whole sky. This sky-integrated flux is the product of the local WIMP num- 
ber density and of the mean speed of the WIMPs with respect to the Earth. It's 
the latter that is modulated. So here we study the time variations of the mean 
WIMP speed with respect to the Earth. (Notice that this is the mean speed, 
not the mean velocity.) 

First we write an expression for the velocity of the Earth. We neglect the 
cllipticity of the Earth orbit and the non-uniform motion of the Sun in right 
ascension (an error of less than 2 days in the position of the modulation maxi- 
mum and minimum). Hence we write the velocity of the Earth in terms of the 
Sun ecliptic longitude X(t) as 

v e (i)=v e [ei sin X(t) - e 2 cos X(t)} , (15) 

where = 27rA.U./yr = 29.8 km/s is the orbital speed of the Earth, and the 
unit vectors ei and &2 are in the direction of the Sun at the spring equinox and 
at the summer solstice, respectively. In galactic coordinates, 

ex = (-0.0670, 0.4927, -0.8676), (16) 
e 2 = (-0.9931, -0.1170, 0.01032). (17) 

The Sun ecliptic longitude X(t) can be expressed as a function of time t in years 
with t = at January 1 as 

X(t) = 360° (t - 0.218). (18) 
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Here 0.218 is the fraction of year before the spring equinox (March 21). 
We give now the mean WIMP speed with respect to the Earth. It is 



(v(t)) = J uf B (u,t)d 3 u, 



(19) 



where /©(u, t) is the WIMP velocity distribution in the rest frame of the Earth. 
In terms of the velocity distribution in the galactic rest frame /(v), 



/©(u,i) = /(u + v +v©(i)). 



So we can also write 



(v(t)) = J | v -v -v e (t)|/(v) d 3 v. 

In the standard halo model, the mean WIMP speed on Earth results 

(v st d(t)) = v (t)g (t), 
where vo(t) — |v Q + v©(i)| and 



.90 (t) = 



+ 1 + 



2x 2 



erf(x ) 



2 1 + zl + xl/3 



x 



with x Q = v (t)/v and z = v csc /v . 

In the SLI halo model with no velocity dispersion, 



<w*)> = £^ 



(20) 



(21) 



(22) 



(23) 



(24) 



with Vi(t) — |vj — v© — v©(i)|. Including a velocity dispersion as described 
above, 

<«gEf(*)) = E^ft(*)> ( 25 ) 



where 



9i{t) 



1 



1 

2x1 



erf(xi) 



(26) 



Xiy/lT 

with Xi = |vj - v© - v©(£)|/t7j. 

We have plotted the mean WIMP speed (v(t)) as a function of time in Fig. 3 
for the standard halo model and the SLI halo model with and without a velocity 
dispersion. The time axis starts January 1 and covers a year. 

Fig. 3 clearly shows that the modulation in Sikivie's late-infall halo model has 
a phase opposite to the modulation in the standard halo model. The maximum 
of the mean WIMP speed on Earth occurs in early June in the standard model 
but in early December in Sikivie's late-infall halo model. With our assumptions 
and approximations, the maximum mean speed occurs June 1 in the standard 
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Figure 3: Annual modulation in the mean WIMP speed on Earth as a function 
of time. The time axis starts January 1 and covers a year. Solid curves in 
SLI halo model (we also include possible velocity dispersions); dashed curve in 
the standard halo model. The phase of the modulation is opposite in the two 
models. 
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model and December 12 in the SLI model; the minimum occurs December 1 
in the standard model and June 11 in the SLI model. We remind the reader 
that the uncertainties in the Sun peculiar velocity and the approximations in 
the velocity of the Earth amount to an uncertainty of a few days in the position 
of the maxima and minima. 

The change of phase in the modulation is independent of the velocity dis- 
persion assumed in the SLI model. A dispersion of 30 km/sec gives very similar 
results to the zero dispersion (see Fig. 3). Even a dispersion of 220 km/s (not 
plotted), similar in magnitude to the average velocity, would give a phase of 
the modulation equal to the zero dispersion case (but the value of the average 
relative velocity would be much larger). 

The local dark halo density in Sikivie's late-infall halo model, obtained by 
summing up the densities in each flow in Table 1, is p= 0.37 GeV/cm 3 . The 
upper bounds on the scattering cross section for a particular WIMP depend on 
the value of the product p(v(t)}. The bounds obtained with Sikivie's late-infall 
model and with a standard model with the same local density would differ by 
the ratio of WIMP velocities in Fig. 3. Otherwise the ratio of the local densities 
in both models should be taken into account. 



4 Recoil-energy spectrum 

The WIMP detection rate per unit detector mass and nucleus recoil energy in 
the range (E, E + dE) can be written as 

*< B >= £f S^"< £ -"- (27) 

Here cto is a normalized WIMP-nucleus cross section; /i — mM / (m + M) is the 
reduced mass of the WIMP-nucleus system (m is the WIMP mass and M is 
the nucleus mass); F(q) is a nuclear form factor function of the nucleus recoil 
momentum q = 2ME; r/(E, t) determines the annual modulation and depends 
only on the WIMP velocity distribution in the rest frame of the Earth /©(u, t), 

T]{E,t) = \ ( dn u f u/ e (u,t)du. (28) 

A J Jy/ME/2/j, 2 

The nuclear form factor F(q) depends on the type of WIMP-nucleus in- 
teractions, namely if they are spin-dependent or spin- independent, and reflects 
the mass and spin distributions inside the nucleus. It may strongly affect the 
counting rates in dark matter detectors, but since it is time-independent it is 
inessential in the analysis of the annual modulation. Therefore in the examples 
below we set F(q) = 1. The recoil-energy spectra we plot should be multiplied 
by the appropriate F(q). 

We now focus on r](E, t) since it is the only time dependent part in the 
recoil-energy spectrum due to WIMP-nucleus collisions. For the standard halo 
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model, 
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(29) 



with w — y / ME/2fi 2 . In Sikivie's late-infall model without velocity dispersion, 
each flow contributes a flat spectrum up to the maximum recoil energy 

£max,> (*) = -gr | Vi - V - V® (t) | . (30) 

The resulting spectrum is a series of steps where each step corresponds to a 
flow, 

VSLi{E,t) = } — — — . (31) 

*f P 2|vi-v -v©(t)| 

Here 9 is the Heaviside function. Notice that the position of the steps given 
by E max ^(t) depends on time, and so the end-point of each step is annually 
modulated. 

In Fig. 4 we plot the recoil-energy spectra at the maximum and minimum of 
the modulation in the mean WIMP speed on Earth. At the lowest recoil-energies 
the phase of the modulation is opposite in the two models; at intermediate 
energies it is the same; at higher energies it is opposite again. To produce 
the figure we have assumed a local halo density p = 0.3GeV/c 2 /cm 3 , a WIMP 
mass m = 60GeV/c 2 , a 73 Ge detector, a WIMP-nucleus cross section a = 
10~ 35 cm 2 , and a nuclear form factor F(q) = 1. Corresponding spectra for 
other detectors can be obtained using the relations E' = (p! //j,) 2 (M/M')E and 

R' = F\q){a' /a ){n/n'fR- 

One may worry that the lowest-energy step does not persist when additional 
flows with smaller velocities are added in the SLI model. We argue that this will 
not happen because while their velocities are smaller and smaller in the galactic 
rest frame they are not smaller and smaller in the Sun or Earth rest frames 
(see column 5 of Table 1). Hence the end-point of the lowest-energy step, which 
corresponds to the flow with the smallest speed relative to the Earth, will not 
change if additional flows above the 20th are added in the model. 

Fig. 5 shows how the recoil-spectrum is modulated at fixed recoil-energies 
E = 15, 25, 35 keV. The change of phase of the modulation is again visible. It is 
also clear that except at the lowest recoil energies the spectral modulation in the 
SLI halo model is no longer well approximated by a sinusoidal function. This is 
due to the modulation of the position of the steps at E maK< i(t) mentioned above. 
If particle flows are introduced in the halo model, a more general analysis of the 
experimental data than that used up to now may be called for. 

The effect of a velocity dispersion in each flow is to blur the end-point of 
each step. The recoil-spectrum in the SLI model with velocity dispersion can 
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Figure 4: Recoil-energy spectrum at the maximum and minimum of the annual 
modulation in the mean WIMP speed on Earth. Step-like curves in the Sikivic's 
late-infall halo model; smooth curves in the standard halo model. The phase 
of the modulation in the two models is opposite at low energies, the same at 
intermediate energies, opposite again at high energies. Notice that the end-point 
energy of each step in the SLI model is annually modulated. 
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Figure 5: Annual modulation of the recoil-energy spectrum at several fixed 
energies. Solid curves in the Sikivie's late-infall halo model; dashed curves in 
the standard halo model. Beyond the end-point energy of the lowest step, the 
modulation of the recoil-energy spectrum in the SLI halo model may be poorly 
approximated by a sinusoidal (cfr. the 25keV curve, e.g.). 
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be obtained using the expression of rj st d(E,t) after taking the limit v csc — ► oo 
and replacing l/A osc -> J^i Pi/P and v eM l v i - v - v ffi (t)|. 

In Sikivie's model there is no virialized halo component. If one should exist, 
as assumed in Ref. 0, the local halo velocity distribution and local WIMP 
fluxes would be a combination of both components. This would result in a 
superposition with some relative amplitude of the fluxes and rates shown in 
Figs. 1, 2, 4, and 5. The phase of the annual modulation would then depend on 
the relative strength of both components. 

5 Conclusions 

The late-infall model j3J of halo formation, modified in recent years by P. Sikivie 
and collaborators ^, |J ^| to include a net angular momentum and axial sym- 
metry around its direction, predict non- virialized flows of dark matter particles 
falling into the galaxy and oscillating in and out many times. Hence it predicts 
a local velocity distribution completely different from the standard truncated 
Gaussian usually assumed. P. Sikivie in Ref. || has given the local velocities 
and densities of the first twenty pairs of flows (the first pair corresponding to 
particles coming for the first time into the galaxy from opposite sides of it, the 
second to those passing for the second time, etc.) in his particular late-infall 
model that has parameters that fit well the halo of our galaxy. In this model 
most WIMPs come from directions in the hemisphere opposite to the Sun mo- 
tion. Namely the average "WIMP wind" velocity on Earth (as pointed out by 
Sikivie himself ||) is reversed with respect to that in the standard halo model. 
Thus, in Sikivie's halo model, the annual modulation of a galactic WIMP signal 
has a phase opposite to the usually assumed. This was illustrated in Fig. 3. 
Moreover (as shown in Fig. 2) many WIMPs in this halo model approach the 
Sun and the Earth from directions above and below the galactic plane, with a 
clear pattern of directions. This would be very important for experiments in 
which directionality could be measured. 

The main message we want to convey is that if a non-virializcd halo compo- 
nent due to the infall of (collisionless) dark matter particles cannot be rejected, 
an annual modulation in a dark matter signal should be looked for by experi- 
menters without fixing the phase a-priori. The detection of an annually modu- 
lated dark matter signal would not only reveal the nature of the dark matter in 
our halo, but the structure of the dark matter halo itself. 
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